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Random matrix models consisting of normal matrices, defined by the sole 

constraint [N T ,N] = 0, will be explored. It is shown that cubic eigenvalue 
repulsion in the complex plane is universal with respect to the probability 
distribution of matrices. The density of eigenvalues, all correlation functions, 
and level spacing statistics are calculated. Normal matrix models offer more 
probability distributions amenable to analytical analysis than complex matrix 
models where only a model with a Gaussian distribution are solvable. The 
statistics of numerically generated eigenvalues from gaussian distributed nor- 
mal matrices are compared to the analytic results obtained and agreement is 
seen. 

PACS numbers: 02.10.Sp, 02.30.Fn, 05.45.+b, 03.65.-w 

I. INTRODUCTION 

Random matrix theory (RMT), has found much success as phenomenological 
models describing a wide variety of physical systems, from discretization of moduli space 
||, to the statistics of the cells in the skin of a cucumber M. RMT is the study of eigen- 
values derived from random ensembles of matrices with stochastic elements specified by a 
probability density, P(M)<iM, in the space of matrices. Most interest is in examining the 
properties of the eigenvalues induced from the transformation to the eigenvalue basis. Early 
work consisted of using real symmetric, hermitian, unitary, and real quarternion matrices 
||. The eigenvalues of which are either real or unimodular, (and termed one dimensional 
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eigenvalues here). 

Recently several groups have begun to consider physical applications of matrix models 
composed of complex matrices. Introduced in the early sixties by Ginibre || it took decades 
for others to consider applications. As the probability distribution is not invariant under 
similarity transformations the diagonalizing parameters must be integrated out by brute 
force. It was found that only a distribution with a gaussian weight (Ginibre Ensemble) 
could be solved. The Ginibre ensemble was later found to exhibit cubic eigenvalue repulsion 
in the complex plane, @ |§. Such models are of interest in characterizations of quantum 
chaos. Cubic quasi-energy level repulsion is a key signature of classical chaos within a 
quantum dissipative system, as defined by Haake 0. 

Normal matrices are discussed in most matrix theory texts, (see ||, 0). Defined by 
the sole constraint that they commute with their adjoint, they have the property of being 
the most general matrix that can be diagonalized by a unitary transformation. The normal 
matrix model was first introduced in showing how the Laughlin wavefunction could be 
modeled by it and offered some generalizations to inhomogeneous fields |TD|. In |TTJ the 
statistics of eigenvalues of random normal matrices were first explored. 

In this paper the universality of cubic eigenvalue repulsion in the complex plane for en- 
sembles of random normal matrices will be shown. In section II the probability distribution 
in the space of matrices, all correlation functions, and the eigenvalue densities will be ob- 
tained. In section III the level spacing statistics of the eigenvalues obtained will be derived. 
The term level spacing is borrowed from the Wigner-Dyson ensembles and will represent 
here the spacing in the complex plane. In section IV eigenvalues of numerical generated 
normal matrices will be studied and compared to the analytical results. 

II. PROBABILITY DISTRIBUTIONS, N-PT. CORRELATION FUNCTIONS 

We begin by defining the joint probability distribution (jpd) within the space of normal 
matrices, 
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Pjv(M) = C exp[-Tr(V{M,M r ))]. (1) 

Only potentials which are hermitian will be considered, V(M, M^) = V(M, Mt)t. As a 
result the joint probability distribution for the eigenvalues will be rotationally symmetric in 
the complex plane. Normal matrices are the most general matrices which can be diagonalized 
by a unitary transformation, M = UZU^, so that the measure and the weight are invariant. 
Thus it is simple to derive the jacobian for the transformation to the eigenvalues. 



Proceeding analogously to the hermitian matrix case [[L2|], the metric in the space of 
normal matrices is defined as, 

(ds) 2 = Tr [dM f dM] 

= Tr(Ud(U t ZU)U t Ud(U t Z*U)U t ) 

= Tr([dS, Z][dS, Z*] + ([dS, Z]dZ* + dZ[dS, Z*]) + |dZ| 2 ) 

= c'£\dS ik (z k -z i )\ 2 + y £\dz i \ 2 1 (2) 

i<k i 

where dS = UdU^ = — dUU^ is antihermitian and c is an overall constant. 
Using the standard Riemannian volume form with this metric we get, 

d/i(M) — ► dfi{z) = (dVL)\/\(z)\ 2 dzf ] A dzf ] A ... A dzf A dz? A ... A dz^ . (3) 

where dfl is the volume of the unitary group U(N), and A(z) is the well known Van der 
Monde determinant, 

N 

iA(*)i 2 =ni%-^i 2 w 

i<j 

If we consider weights which are invariant under unitary transformations we can factor 
dfl out, 

N 

P N {z) = C'\ A(^)| 2 exp[- £ V(zi, z*)}, (5) 

i=i 

where C is a normalization constant. 

Considering a gaussian weight, V = \z\ 2 , we notice that the jpd is identical to Ginibre's 
ensemble of complex matrices |§ . The model with this weight function is termed the Normal 



Gaussian Ensemble (NGE). It is clear that gaussian ensembles of complex matrices and the 
NGE are identical. Any result derived here for the NGE will also be valid for Ginibre's 
ensemble of matrices. This is of little surprise since TrfM^M] = Tr[MMJ]. 

The difference occurs when considering other ensembles. For complex matrices the pa- 
rameter /3, (defined as the power of the Van der Monde determinant appearing in the jpd), 
is a function of the weight. In fact, only models involving gaussian weights have ever been 
studied as this is the only nontrivial model which allows a separation of radial and angular 
parameters. For normal matrices (3 = 2 regardless of the weight and it is possible to study 
a wide variety of ensembles. 

For matrix models with one dimensional eigenvalues the traditional method of analysis 
employs a basis of orthogonal polynomials. These allow the reduction of the determinant 
upon integration of a number of eigenvalues || . For complex eigenvalues we can also intro- 
duce an orthogonal basis. If a polar basis is chosen, (and Pn(z) is rotationally symmetric), 
it is simple to verify that the z's themselves form an orthogonal basis. A basis of orthogonal 
monomials, {</>}, is defined, 



= T"> (6) 

where Ni = J °° dr 2 r 21 exp[— V(r)], is a normalization constant. 
The jpd can then be expressed as, 

P N {z) = C'e-^=i v ^detK N {z hZj )\ i>j=1 , N , (7) 

where the "kernel" is, 

JV-l 

K N (zi,Zj) = M^tifa)- ( 8 ) 

1=0 

The weight function is deliberately factored out of the kernel since all of the angular depen- 
dence is contained within K. 

Similarly to hermitian matrix models, when N — n eigenvalues are integrated out of the 



jpd, the determinant shrinks by N — n columns and rows, [IT]. This property allows the 
n-point correlation function to be derived easily, 



i\n r l 
i^fci, • • • , z n ) = — — / J] -de i dr 2 i P N (z), 

\ /" i=n+l 



= e -Er =1 nn) defjft:Ar ^.^.)|. i=1)R . (9) 

The prefactor Tj^z^j is due to the ordering of eigenvalues. The 2 point correlation function 
and eigenvalue density are, 

e -V(z)-V(w) 

R 2 (z,w) = \p(z)p(w) - K(z,w)K(w,z)), 

N-l I 12 

P (z) = e-^W £ ^, (10) 

In deriving expressions for the correlation functions orthogonal monomials provide a more 
economical method. However, there are no simple recursion relations among the monomials 
as there are for the orthogonal polynomials. No recursion relation implies no Christoffel- 
Darboux formula and no asymptotic analysis. However, we can find analytical results for 
finite N. It is not necessary to examine the asymptotic form of the monomials to get large 
N forms. 

In table | various ensembles are defined and their corresponding eigenvalue density and 
two point correlation function are given. It should be noticed that several of these ensembles 
allow global closed forms in the limit iV — > oo. 



III. LEVEL SPACING STATISTICS 

In order to study the properties of random matrices the statistical properties of the 
spacings between eigenvalues, (or level spacings), are examined. These statistics can be 
used as a definition of quantum chaos. The distinction of classically integrable and chaotic 
systems can be seen in a semiclassically quantized system as the transition from Poisson 



to Wigner-Dyson statistics of the nearest neighbor distribution of the energy levels It 
has been shown J7[ that if there is dissipation in the system the level spacings undergo a 
similar transition from a Plane Poisson distribution, (random points in a plane), to Ginibre's 
distribution generated by an ensemble of random complex matrices. 



The distributions for eigenvalues in the complex plane will require the definition of sta- 
tistical quantities analagous to those defined for Hermitian matrix models 0. 

• The probability that no eigenvalue lies within a circle of radius s centered upon the 
point w, the gap distribution, is denoted as E(s,w). A related quantity E(s,w,e) is 
defined as the probability that no eigenvalues lie within an annulus of outer radius s, 
inner radius e, and centered on the point w. 

• The probability that no eigenvalues lie within this same circle but one or more eigen- 
values lying on the perimeter is denoted as F(s, w). 

• For an eigenvalue lying at the center of the circle and no eigenvalues within the circle 
nor on the edge, this probability is denoted by H (s, w). 

• The level spacing probability, p(s,w), is the most quoted quantity and is defined as 
the probability that one eigenvalue lies at the center and one ormore on the perimeter. 

By taking differential areas at the center and at the perimeter of this circle, we can obtain 
all of the probabilities from the gap distribution. The argument is analagous to Mehta's for 



one dimensional eigenvalues and can be found in |Tl 

N 

E(s,w)= l{dz}P N ({z})l[[l-x(z i: s : w)], (11) 

i=l 

F(s,w) = ~E(s,w) (12) 

H(s,w) = --^E(s,w,e)\ £=0 (13) 

p(s,w) = -— H(s,w) = (14) 

--r^F(s,w,e)\ e=0 = —-^E(s,w,e)\ e=0 . 
ae* as ae 

Where, \ is the characteristic function and e is taken as the radius of an infinitesmal area 
centered at the center of the circle. These relations allow an elegant derivation of the level 
spacing distribution. 
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All of the probability distributions considered here are symmetric in the Zi, hence the 
derivation of the gap distribution will be the same as in ||, (see also JO]]). For the circle 
centered at the origin we have the following, 



e( S ,o)= n 

i=0 

where x = \z\ 2 . 



/ °° dxexp(— V(x))xi 



(15) 



As a check on our results, we find the spacing distribution for the NGE and compare 
with previously known results. In this case the probability distribution is translationally 
invariant, hence the result flT5|) is valid for all w. The result for small s is, 

p(s)^2s 3 -s 5 + ±s 7 -^s 9 --- (16) 

This result is the same as for complex matrices derived by Haake J7|. The level spacing 
statistics at the origin for other ensembles of normal matrices can be found in [11]. Here we 
are concerned with obtaining a universal result for the spacing distribution. 

For ensembles other than Gaussian, the distributions are not translationally invariant and 
hence it will be necessary to resort to asymptotic approximations. Shifting the eigenvalue 
coordinate, z = z — w, does not affect the measure or Van der Monde determinant. The gap 
distribution can then be expressed as, 

E(s,w) = J] / drldO^l - X i(s 2 )]P N (z), 

N N 

= / K^}{1-EX^ 2 )+ E Xi(s 2 ) Xj (s 2 ) -■■■}Pn(z), 



i=l i^3=l 
N „ N 



= 1 - E / d 2 ZiXip{zi) + E / d 2 z i d 2 z j XiXjR2(z i , Zj) , (17) 

i i<j 

where f(z) = f(z — w) and Xi is the characteristic function for the i th eigenvalue. In the 
limit s — > 0, and keeping only the first two terms we obtain, 

]im E(s,w) w 1 - Ntts 2 p(w). (18) 

The level spacing distribution, p(s,w), can be found via equation (|I3]) . First the gap 
distribution E(s,w,e) needs to be derived. The derivation is the same as before except for 
a modification to the characteristic function. Starting from (|L7|) we have, 



N . N . 

E(s, w,e) = l- ^2 J d 2 ZiXip(zi) + Y. J d 2 z i d 2 z j XiXjR2(z i , zf) , (19) 



N „ N 

E 

where now x% — $(s 2 ~ r i)@{ r i ~ g2 )- As an intermediate step we find, 



d 

d7 2 

d -- a 



H(s,w) = -—E(s,w,e)\ e=0 , 



f dr 2 [ d8p(z + w) 
de l Je 2 J 

N 2 + N d r* 2 r° 2 drfdr 2 r r Jn Jn _ . 

J 2 J 2 47 j de 1 de 2 R 2 {z l + w 1 z 2 + w) + ---\ e=0l 



2 de 2 Je* 7,2 

< 2 dr 2 



r civ r 

= -Ntx P {w) + n(N 2 + N)J o -± J de 2 R 2 (w, z 2 + w) . (20) 

The last term exploits the symmetry in R 2 . For the level spacing we have, 

p(s, w) = -n(N 2 + N)s J d6 2 R 2 (w, se id2 + w) ■ ■ ■ (21) 

In order to find the lowest order contribution to p(s,w) the lowest order of R 2 will need to 
be found. For a model with weight e~ v ^ we have, 

R 7 , „a _ -v {w) -v (z+w) f f r 2 J(r 2 + zw* + wz* + r 2 w r _ {zw * + r 2J {w z * + ^jm 

1=1 m=0 J^l^m iV^V m 

The normalization factor, Ni is model dependent. It is easily verified that terms with I = 
and I = m vanish. All higher order functions will have the same property. 



The remaining angular integral over d6 2 in (^l|) insures that no terms linear in z or z* 
appear. It is also apparent that the leading term in r w cancels, (terms to order z°). This 
yields the lowest order of R 2 to be s 2 . Combined with equation (|14"D , it is found that the 
nearest level spacing is cubic in s as s —> 0. A few exceptions occur when the weight function 
vanishes as z vanishes. Specifically the Jacobi, Laguerre, (a > 1) and generalized Gaussian 
(a > 1) will have a higher order repulsion. Our final result is thus: 

Normal matrix models with normalizable hermitian weight functions have eigenvalues 
experiencing a minimum of cubic repulsion in the complex plane. 

It is also worthwhile to examine the case of real normal matrices. Here the constraint 
being that the matrix commutes with its transpose. Such a matrix is diagonalizable by a 



S 



unitary transformation and will in general have complex eigenvalues. Retracing the steps 
above it is easy to see that the eigenvalues will also exhibit a minimum of cubic repulsion 
in the complex plane. 



IV. NUMERICAL RESULTS 



Whenever discussing random matrix theory it is always of interest to include numerical 
simulation of large random matrices. These are used to verify the analytical properties and 
supply support for asymptotic forms, conjectures, etc. When dealing with normal matrices 
however, we are in a predicament. Consider the following problem: 

• Write down an arbitrary 2x2 normal matrix that is not Hermitian, Unitary, nor Real. 

It is not so simple. If one was asked to write down a hermitian matrix, it could be done 
without thinking twice. For a normal matrix there is no simple parameterization. It would 
most likely require the commutator to be calculated and a solution found by trial and 
error. Likewise, if you want a computer to generate normal matrices you must give it a 
parameterization. For a 2 x 2 you can parameterize it as follows, given that n n * ^ n 2 2*, 
(remember there are N 2 + N = 6 independent real elements within a normal matrix), 



For arbitrary N the constraint equations become nonlinear and highly coupled. We will 
choose the following parameterization of normal matrices, define the independent elements 
as nij, (those along the diagonal and above the diagonal). Define the elements that are 
constrained as fji, (those which lie below the diagonal). For the ji th element the constraint 
equation is, 




(23) 



fji 





i-l N 



im n jm 



1 



i-l N 




(24) 
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Notice that the term in the first bracket contains only independent elements. With this 
parameterization the probability of obtaining normal matrices which are hermitian, unitary, 
orthogonal, or real symmetric matrices becomes infinitesimal for large N. The subspaces for 
these types of matrices within the space of normal matrices are of much lower dimensionality. 
There are many solutions to this nonlinear equation. Another method must be devised to 
obtain a normal matrix. It is also found that due to the nonlinearity, only normal gaussian 
matrices can be generated. (Others not respecting the hermiticity and normalizibility could 
be analyzed but is not of interest here). 

Instead of directly attempting to produce random normal matrices we begin with a 
random complex matrix and use an optimization algorithm to obtain a normal matrix. The 
following algorithm is the simplest to implement but by no means the most efficient. 

1. Parameterize the normal matrix as before: 

Elements along the diagonal and above the diagonal are independent parameters. 
Elements below the diagonal will be optimized. 

2. Using the Box-Muller technique, generate random gaussian variables and construct a 
200 x 200 complex gaussian matrix. This is the initial condition. 

3. Calculate the commutator and define the following: error function, tr[N,N f ] 2 = E. 
The error function will be minimized below a choosen threshold, (a value of 0.2 is 
used) . 

4. Use gradient descent to adjust the constrained elements. If the value of E is below 
the choosen threshold then pass the matrix to the diagonalization routine. If E hasn't 
converged after 500 iterations, then reset the matrix and go to step 2. 

The results of a preliminary numerical study based on this algorithm are now presented. 



A full discussion of the algorithm used is to be found elsewhere, (TT|. Due to limited 
computational resources and time, only 50 matrices of size 200 x 200 have been analyzed. 
The level spacing statistics for the generated matrices is compared to the analytical result. 
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Before attempting to generate large normal matrices the algorithm was checked with 
2x2 matrices. As these were easily produced, 11,000 were generated and compared to the 
exact analytical result, 

,2 



y K ' 10.4683 V 4 ; 



(25) 



The eigenvalues are unfolded so that the mean spacing is unity, this is achieved by 
rescaling the spacing and renormalizing the spacing distribution so that the mean spacing 
between eigenvalues is unity, 



S = 1.6622s, 
p(S) = p(s)/10.4883. 

Plots of the analytical form and the unfolded eigenvalues are displayed in figure |l|. 



(26) 




1 1.5 

Spacing, S 

FIG. 1. The level spacing distributions, p(S), for N = 2, analytical form and 11,000 numerically 
generated normal matrices. 

For normal matrices of size 200 x 200 the exact analytical form for the spacing distribution 
can easily be found from equations (|T]])-(|T4"|), 



199 e 1 (s 2 ) 
p 2 oo( S ) = (199-E 9 1 j 



199 



)2sY[[e- s2 ei (s 2 ) 



(27) 



i=l 
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where e q (s 2 ) = J2i=o if- Under the optimization algorithm it is found that many eigen- 
values "leak" out into the complex plane. Adding a further constraint to keep the eigenvalues 
within a radius of y/2aN — 1 « 20 lengthens the generation time dramatically. Therefore, 
for this preliminary study the eigenvalues were "pruned", those falling outside r = 18.1 
were neglected, (the number was choosen so that edge effects would be eliminated as well). 
Overall, about one third of the eigenvalues were pruned. The results are displayed in figure 

! 



P(S) for N - 200, 50 matrices 




Spacing, S 

FIG. 2. The level spacing distributions, p(S), for N = 200, analytical result (with eigenvalues 
unfolded), and 50 numerically generated 200 x 200 normal matrices. 

The slight deviation of the numerical and analytical results for large S is due to the 
leakage of eigenvalues. The small S behavior is not affected by the leakage. It is found that 
changing the initial conditions has little effect on p(S). It is also observed that lowering the 
the threshold for the error function by a factor of ten has a negligible effect. 

V. CONCLUSION 

It has been shown that cubic repulsion of eigenvalues obtained from random normal 
matrices experience a minimum of cubic repulsion in the complex plane. Even though normal 
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matrices are subsets of complex matrices we have demonstrated that Gaussian normal matrix 
models yield the same results as the Ginibre ensemble of complex matrices. One can define 
ensembles of complex matrices with weights invariant under similiarity transformations, 
(see appendix B of JIT]]). These ensembles will fall under the /3 = 4 and experience a 
minimum of quartic repulsion. (Getting analytical results for the statistics of the eigenvalues 
is difficult however). Ensembles of complex matrices with weight functions invariant under 
unitary transformations other than the Gaussian are not factorable and are not amenable to 
analytical analysis. For normal matrices we can define identical ensembles to complex models 
invariant under similiarity transformations and as well define a wide range of ensembles 
invariant under unitary transformations, (as demonstrated in this paper). The special case 
of real normal matrices was also shown to experience cubic repulsion in the complex plane. 
As well, there are even more significant differences between random normal matrices and, 
more traditional, random matrix models composed of hermitian, unitary, and real symmetric 
matrices, themselves being subsets of normal matrices. Although numerical generation of 
normal matrices is not as straightforward as for other types of matrices we have shown that 
analytical results are often much easier to obtain. What numerical information we could 
obtain does agree with the theoretical result. 



It is a little disappointing that universality, in the sense of Brezin and Zee ||14|| , is not 
found for normal matrix models. The simple scaling arguments that lead to universal kernels 
for other types of matrices do not work here. It may be possible to rescale S to obtain a 
universal spacing distribution P(S). This allows the potential of analytically exploring the 
statistics of complex eigenvalues which was not possible before. Another interesting feature 
is the closed form expressions for the two point correlation functions in NGE, Laguerre 
(a = 0), and Legendre ensembles in the infinite matrix limit, (see Table |). 

In terms of physical applications it should be noted that results using the Ginibre en- 
semble of complex matrices can be duplicated with the normal Gaussian ensemble. Such 
work includes, modeling cellular structures by Le Caer and Ho H , two dimensional quantum 



gravity by Morris and Haake's work on generators of quantum dynamics for systems 
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with dissipation 0. It has recently been shown, |L£]], that Haake's generator of dissipation, 
D, in [IT] |7|, is properly represented by a normal matrix rather than a complex one. This 
result strengthens the idea that the eigenvalues of these generators universally exhibit cubic 
repulsion when dissipative dynamics is present. 

The results obtained in this work give a more complete picture of the statistics of eigenval- 
ues from random matrices. The normal matrix model, a link between the Dyson ensembles 
and the Ginibre ensemble, has a rich structure where interesting results can still be found. 
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TABLE I. Ensembles of random normal matrices defined by their probability distribution. Also 
shown are the corresponding eigenvalue density and two point function. For the NGE, Laguerre 
a = 0, and Legendre ensembles the closed form large N result has been shown. 
Ensemble Probability Distribution Density, p p(ri)p(r2) — R2 



NGE g-TrM'M 1 J^exp[-|z-™| 2 ] 

NQE e -T,(MtM- 9 MtMtMM) e_f_f^_ ^N-l _r^_ ^f^ 4 K N { Zl , Z 2 )K N {z 2 , Zl ) 

Tfl „ llprrp _-Tr((MtM)3+Sln(MtM)) e -y°+i r 21 v iV-l (n^^e-n-a^*^*^™ 

I^agueiie e 2?r 2^l=0 r[2i+a+2| ^1=0 4n 2 r[2l+a+2]T[2l+a+2] 

Laguerre, a = e -Tr(Mt M )i) |_L sinh(r ) sin h( v ^4) sinh(yi*^) 

Legendre 1, (evals G unit circle) ^rhpiYi rr=^pl 0- ~ z \ z % 



Tacnbi r/p/M - nVrtM^la (l-r)° ^N-l r 21 1 Zl l z* l z* m z^ 

JdGOUl UCL[1 (W± iv±j j 2w 2^1=0 B[a+l,2i+2] 4^^1=0 B[a+l,2l+2]B[a+l,2m+2] 

2 2 

fVr, Pnnaainn p -Tr (M*M-% ln(MtM)) e^V V AT-1 r 2 ' v iV-l (riT- 2 ) n e - r l r 2 Zl l z* l z* m z^ 

uen. uaussian e ^ 2^z=o r[z+f+ij ^«=o ^r[/+f +i]r[ m +§+i] 
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